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Abstract 

Strong correlation effects in classical and quantum plasmas are discussed. In particular, Coulomb 
(Wigner) crystallization phenomena are reviewed focusing on one-component non-neutral plasmas 
in traps and on macroscopic two-component neutral plasmas. The conditions for crystal formation 
in terms of critical values of the coupling parameters and the distance fluctuations and the phase 
diagram of Coulomb crystals are discussed. 
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I. INTRODUCTION 



Coulomb crystals (CC) - a periodic arrangement of charged particles - are omnipresent 
in nature, from astrophysical systems (interior of dwarf stars, Refs. to laboratory 

systems (trapped ions, see e^. Refs. 0, 4|, plasmas in storage rings, e.g. Refs. [s], ^1 or 
dusty plasmas, Refs. [3, Isl, IsJ, to name a few examples, for an overview see Ref. jlO|). 
CC add an interesting new species to the large family of crystals in condensed matter, 
chemistry or biology, for an (incomplete) list, see table [B We will distinguish CC from 
"traditional" crystals (including molecular or ion crystals or metals) by 1. the governing 
role of the Coulomb interaction (in contrast to crystals of neutral particles) and 2. by 
the elementary character of the constituents (in contrast e.g. to the complex ions forming 
the lattice of a metal) These two properties bring the CC into the area of plasma 

physics rather than condensed matter physics, because it is the strength and long range of 
the Coulomb interaction which dominates the many-particle behavior in these systems, the 
crystal symmetry, stability and melting properties. 

The research on CC originates in solid state physics. More than seven decades ago Wigner 
predicted, using the jellium model, that electrons in metals would form, at low density, a 
bcc lattice, see Ref. [l2|. A second line of research grew out of the field of classical strongly 
coupled plasmas. There it was predicted, by computer simulations, that a one-component 
Coulomb or Yukawa model plasma (OCP) in three and two dimensions would crystallize 
at sufficiently high density and/or low temperature, e.g. Ref. [10|]. 3D Coulomb crystals 
show a bcc symmetry whereas Yukawa crystals have a bcc and a fee phase, Ref. 13|. In 
contrast, the ground state of 2D crystals has hexagonal symmetry. However, jellium and 
OCP are models assuming that the charge species forming the crystal coexists with a second 
neutralizing one which forms a static homogeneous background which does not influence the 
crystal. Such systems do not exist in nature. In real two-component plasmas crystallization 
is very different. One important effect is weakening of the Coulomb interaction by dynamic 
screening. Moreover, the attractive force between different species will favor recombination, 
i.e. formation of bound states. This will, obviously, strongly reduce the Coulomb coupling 
and may even prevent crystal formation. Nevertheless, CC formation in a two-component 
plasma (item AIII. (c) in table [T]) is possible and will be discussed below in Sec. IVII 

But before that we consider the second possibility to achieve Coulomb crystallization: 
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one-component (non-neutral) plasmas which are stabilized by an external "trap", such as 
an electric potential, cf. item BII in the table. This principle has been successfully used in 



HQ 



for an overview see Refs. 



19|, 



experiments with ion crystals, e.g. Refs. j^, |4|] and dusty plasmas, e.g. Refs. 



14 



15 



20|, and is expected to function also with electrons 



in semiconductor quantum dots, Ref. [21|. Naturally, the existence of the trap may have a 
strong influence on the crystal properties. For example, a spherically symmetric trap will 
favor crystals forming concentric rings (in 2D) or shells (in 3D). This gives rise to interesting 
symmetry effects, including magic (closed shell) configurations, e.g. Refs. |22|, |23|, l2J, [25 1 



familiar from atoms and nuclei and coexistence of shells and bulk behavior in larger systems, 



Ref. 



26|. 



Coulomb crystals may not only consist of classical "point particles" but also of quantum 
particles which have a finite extension (electrons in quantum dots, ions in compact stars 
etc.) which is of relevance for the properties of CC and is crucial for the phase diagram. 
Since the issue of quantum plasmas has come into the focus of recent research again in the 
context of laser plasmas 27|] and astrophysics 281] we will consider the influence of quantum 
effects in some detail. In this paper we study some general properties of Coulomb crystals. 
Starting from the theoretical description, in Sec. [TTl we continue with two typical examples 
of classical and quantum crystals in traps (Sees. Illlt lIVp . This is followed by an analysis of 
the melting point. Sec. |Vl after which the special situation of CC in neutral plasmas (Sec. 
VB is discussed. 



II. MODEL AND PARAMETERS 



The Hamiltonian of a system of particles with mass rrij and charge Cj interacting via a 
statically screened Coulomb (a Yukawa) potential is given by 



H 



N 

E 

i=l 



N 



2m. 



^2 ^ J ^ 



(1) 



where r^j = |ri — rj| and e denotes a static background dielectric constant, which is of the 
order of 10 in case of an electron-hole plasma in a semiconductor; in a plasma, e = 1. The 
case of a pure Coulomb system follows in the limit of zero screening, k — > 0. In the case of 
trapped systems a confinement potential V{r) is included which will be assumed isotropic 
and parabolic, i.e. V{r) = ma;^r^/2. The limit of an unconfined system is achieved by 



TABLE I: Coulomb crystals (CC) in the world of crystals (incomplete list). CC variants are 
A.III.b, A.III. (c), B.II. (a) and B.II. (b). ICS (2CS) denotes one (two) component systems, OCP 
- the one-component plasma model containing ions plus a homogeneous static neutralizing electron 
background. 



A. Unconfined (macroscopic) crystals 



B. Confined crystals (1, 2 or 3D traps) 



/. ICS with attractive interactions 
neutral particles (e.g. Lennard-Jones, 
Morse potentials) 

"normal" solids, rare gas clusters etc. 



/. ICS with attractive interactions 
confinement not necessary, see A I. 



//. ICS with repulsive interactions 

a. transient "Coulomb exploding" crystals 

b. charges on surfaces of finite systems 
(e.g. electrons on helium droplets) 



//. ICS with repulsive interactions 

a. classical: ions, dust particles 

b. quantum: electrons in quantum dots 



///. Periodic confinement 

e.g. particles in optical lattices 



///. 2CS 

a. "normal" crystals: ionic crystals, metals etc. 

b. OCP model (ion Coulomb or Yukawa crystal) electrons in bilayers, superlattices etc. 

c. TCP crystals (electrons, nuclei, holes, positrons) 



letting uj ^ Q. In thermodynamic equilibrium the system properties are determined by the 
canonical probability distribution P or, in the quantum case, by the density operator p 



P{E) = ie-^^ 



(2) 



where (3 = l/Zc^T is the inverse temperature, and E denotes the total potential energy, i.e. 
the 2nd plus 3rd term of Eq. ([T]) . 

Despite their different form of appearance, all Coulomb (Yukawa) systems exhibit similar 
fundamental properties governed by the strength of the Coulomb (Yukawa) interaction which 
is measured by dimensionless control parameters: the coupling parameters Ta, Vga and \a 



of particle species "a" and the quantum degeneracy parameter Yn,. These parameters are 
determined by the ratio of characteristic energy and length scales 29|, |30[ |: 



• Length scales: (i) f - average inter-particle distance, f ~ n~^^^ [n and d = 
1,2,3 denote the density and dimensionality of the system, respectively), (ii) A - 
quantum-mechanical extension of the particles. For free particles we have A^^*^*^ = 
h/ \/2TimakBTa (DeBroglie wavelength), for bound particles A is given by the exten- 
sion of the ground state wave function, A^"""*^ = 2TTaB. (iii) as - relevant Bohr radius 
Or = — where m~} = + m7^. (iv) an„ - effective Bohr radius of an OCP: 
a = ' 

• Energy scales: (i) {K) - mean kinetic energy, which in a classical system is given 
by {Ka)ci = ^ksTa, whereas in a highly degenerate Fermi system {Ka)qm = f-^Fa 
holds [Ep = ^^(37r^?2)^/^/2m denotes the Fermi energy]; (ii) (f/^^) - mean Coulomb 
energy, given for free and bound particles by {Uf')f = and {U^^)b = 1^2^^ = 
Er (Rydberg), respectively. Analogously, the mean Yukawa interaction energy is 
estimated by (Uy) f = e~'^'^ {Uc)f- 

• Dimensionless control parameters: The quantum degeneracy parameter Xa = ''^aAf ~ 
{Aa/faY divides many-body systems into classical (x < 1) and quantum mechanical 
ones (x > !)• The Coulomb coupling parameter is the ratio \{Uc)\/ (K) . For classical 
systems Fq = \{U^"-)\/kBTa results, whereas for quantum systems the role of F^ is 
taken over by the Brueckner parameter, = fa/asa ~ \ {Uc"')\/ Epa- The relation to 
the parameter, = f/aB, familiar from atomic units is = r^afl — ^ ] . Similarly 
one can introduce coupling parameters for Yukawa systems and of different species. 

In a two- component plasma different masses and charges of the species may give rise to 
unequal coupling and quantum degeneracy of the species. In particular, in a dense electron- 
ion plasma classical ions and quantum electrons may coexist. Analogously ions maybe 
strongly coupled while the electrons are only weakly coupled, see Sec. IVII The ratio of the 
degeneracy parameters scales as Xa/Xb = {'^nb/i^aY^'^ 1 whereas the ratios of the coupling 
parameters are given by Fa/Ft = {ca/ ChY^^^'^ and Vsa/rsb = {nia/mb){ea/ebY^^^'^, where 
local charge neutrality, UaCa = UhCb, has been assumed. 



5 



III. CLASSICAL COULOMB AND YUKAWA CRYSTALS IN TRAPS 





4 




I ' 1 ' 1 ' 1 
K=0 




K=l 
































'v. 






-> • V 1^. .1- 
', ') i; .". 

" . Ji •* 







' ^ 




i ^/ ? 



1 2 01234012345 

radius [mm] r/r 



FIG. 1: Radial particle distribution for N = 190 particles given in cylindrical coordinates. Left: 
experiment, right two figures: simulation results with Coulomb [k = 0), and Yukawa {k = 1) 



potential. The length unit in the right two figures is given by Eq. (Uj), from Ref. 



24| 



Coulomb crystallization in a spherical 3D geometry was first observed for ultra-cold ions 
in Penning or Paul traps [sl. A second candidate are ions created by ionization of cooled 
trapped atoms. Recent simulations 3l[ show that the expanding ions might crystallize 
if they are properly laser cooled durin g th e expansion. Finally, so-called "Yukawa balls" 



have been observed in dusty plasmas l3, 12^, see Fig. [T] Their theoretical description 
is again based on the Hamiltonian ([T]) (for an overview on earlier theoretical results and 
simulations see [l^). This model has, in fact been shown to correctly describe the dusty 
plasma measurements {24, 33|: 3D concentric shells with the populations A^^ being sensitive 
to the screening strength k. With increasing k the reduced repulsion leads to an increased 
population of the inner shells, cf. the table. The quality of the experiments is so high 
that the shell populations can be measured accurately, allowing for comparisons with the 
simulations. In fact, very good agreement is found for htq 0.6, cf. Fig. [21 which shows 
the relevance of screening effects in these confined dusty plasma crystals. Furthermore, 
screening has an important effect on the average radial density profile of these crystals. In 
contrast to Coulomb crystals, where the density is approximately const aiit, with increasing 
K there is an increasingly rapid decay of the density towards the surface 3^, l35 |. 

As in the 2D case closed shell configurations and a "Mendeleyev table" exist (see, e.g.. 



Refs. 



23 



24 



361]). The dependence of the crystal stability on the number of particles can 



be seen from their melting temperatures. For example, the closure of the first spherical 
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shell occurs at = 12, which gives rise to a particularly high crystal stability (high melting 
temperature), cf. Fig. HI 
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FIG. 2: Number of particles on the shells of Yukawa balls with different N and k. Table 
contains experimental (last line) and theoretical shell configurations for N = 190. A'^i . . . denote 
the particle numbers on the i-th shell beginning in the center. Figure shows the shell populations 
for 40 experimentally observed Yukawa balls (symbols) and molecular dynamics simulation results 



for several k values 
units of Eq = e^/ro. 



24( 1 ■ K is given in units of Tq defined by muj Tq = e /rg, temperature is in 



IV. QUANTUM COULOMB CRYSTALS IN TRAPS 

When the trapped CC is cooled, eventually the DeBroglie wavelength A will exceed the 
interparticle distance and quantum effects will become relevant. While for ion crystals this 
may require sub-microkelvin temperatures this regime is easily accessible with (the much 
lighter) electrons in nanostructures. At the same time, there quantum crystal formation 
and detection is hampered by impurities and defects. Therefore, the results shown below 
are obtained by means of computer simulations. The density operator ([2]) with the 2D 
Hamiltonian ([1]) is evaluated by performing first-principle path integral Monte Carlo (PIMC) 



simulations, for details see Refs. 



37 



38| . Results for the probability density of 19 electrons 
in a 2D harmonic trap are shown in Fig. [31 We observe a shell structure similar as in 
the classical case. However, the particles are now not point-like but have a finite extension 
and an elliptic shape which minimizes the total energy. When the system is compressed by 



increasing c^;, the wave functions of the electrons start to overlap - first within each shell, 
cf. central part of Fig. [3l and finally also particles on different shells overlap giving rise to 
a quantum liquid state. This process of quantum melting occurs even at zero temperature, 



giving rise to an interesting phase diagram of quantum CC 



2l|, see also Sec. El 
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FIG. 3: (Color) 19-electron quantum Wigner "crystal" (left), radially ordered crystal (center) and 
mesoscopic fermionic liquid (right). From left to right quantum melting at constant temperature 
occurs. Dots correspond to the probability density p of the electrons in the 2D plane which varies 
between pmax (pink) and zero (red). 



V. CONDITIONS FOR CRYSTAL FORMATION 

Phase coexistence is determined by the equality of the thermodynamic potentials (such 
as the free energy) in the two phases which often requires very accurate and expensive 
calculations. In macroscopic plasmas there exist many alternative criteria for crystallization: 
peaks of the specific heat, sufficiently strong modulations of the pair distribution or the static 
structure factor and so on. These quantities yield practically the same melting point, for 
an analysis see e.g. Ref. 391]. In contrast, in trapped systems, in particular, when the 
particle number is reduced, the results for the phase boundaries may strongly depend upon 
the chosen quantity and the way it is computed. It turns out that, for the class of systems 
described by Eq. ([!]) two quantities are particularly useful to localize the melting point: 
critical values of the coupling parameter and of the distance fluctuations of the particles 
around their equilibrium positions. We mention that, for very small systems, recently a more 
appropriate quantity has been proposed - the variance of the block-averaged interparticle 



distance fluctuations, see Ref. [40|]. 



8 



A. Critical values of the coupling parameter 



Let us start with the simplest case of Eq. ([I]) - a macroscopic classical plasma {u = 0) 
containing a single charge component. We can rewrite the ratio of energy and temperature 
which determines the thermodynamic properties, cf. ([2]), as 

N —Rf- 

(3E = /(r, = r V ^— ^ with fij = ^, R = Kf. (3) 

l<j<i 

For Coulomb systems (k = 0) f3E is characterized by a single parameter, the coupling 
parameter F, i.e. different Coulomb systems (containing different types of particles, hav- 
ing different temperature or density) are expected to show the same behavior if they 
have the same values of F. In particular, as was revealed by simulations, CC occurs 
at Fcr- ^ 175 in 3D and F^^ ^ 137 in 2D. In a Yukawa OCP {k > 0) the effect of 
screening suggests to introduce Fy(fi;) Fe"'"", however, this does not correctly repro- 
duce the K— dependence of the melting curve. The reason is that melting is not deter- 
mined by the absolute value of the energy but by the energy contribution of particle 
fluctuations around their ground state positions rjo- Expanding ([3]) around rjo, defin- 
ing = Tij — fijo and taking into account that the first derivatives vanish we obtain 
(3AE = (3{E -Eo- E,,m) = T Ef<, ^ (l + ^^^J0 + ^) e"^^"-" + . . . . The dots denote 
terms with mixed derivatives and higher order terms, and Eq and Ecom are the energy in 
the ground state and of center of mass excitations (which are not relevant for the melting), 
respectively. 

For the case of two particles, this expression can be written in a Coulomb-like form, 
PfgAE = Fy(K).^^, with the Yukawa coupling parameter Fy(/€) = Fe~'^''"[l + Kf + (Kf)^/2]. 
Assuming that, at the melting point, the critical coupling parameter is universal (3D case), 
Fycr = 175, the phase boundary of the (bcc) crystal in the T — n plane is approximated by 
Fcr(fi;) = 175 ■ e'^'^[l + nf + [Kf)"^ /2]^^. Interestingly, simulations have shown that this result 



no 



ds reasonably well not just for small particle numbers but also in a macroscopic system 



Consider now a classical crystal in a trap. Here the density is externally controlled by the 
trap frequency u which determines the mean interparticle distance. The basic properties 
are best illustrated for two particles. The ground state is obtained from the minimum of 
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the relative potential energy (jl]) with the result 







1 + Kro 



m 2 ~ '"Oc- 



(4) 



Eq- Q yields the two-particle distance, tq^k), as a function of the distance in an unscreened 



system, roc 



24|]. In analogy to the macroscopic case we introduce a Coulomb coupling pa- 



rameter, r2 = / {kBTro). The corresponding coupling parameter for Yukawa interaction, 
T2Y1 again follows from expansion of the energy around the ground state, (Sr^AE = r2y(/t)^^ 



with the result 



2J] = f^^o _^ f^^^ _|_ slightly differing from the above expres- 



sion. In a similar way, the ground state and effective coupling parameter can be defined for 



any particle number, but this has to be done numerically 



23 



42l |. The results are strongly 



A^— dependent due to the importance of shell filling and finite size effects. This leads to 



strong variations of the crystal sta^ 
left part of Fig. m see. e.g. Refs. 



jility with as can be seen in the melting temperatures. 



43 



4|. 



Consider now a macroscopic quantum OCP. We rewrite the Hamiltonian ([T]) in dimen- 
sionless units 



H 



N 



2E 



R 



(5) 



which depends on the quantum coupling parameter and temperature separately, leading to 
a more complex behavior than in a classical OCP where only one parameter F exists. The 
existence of three energy scales - quantum kinetic energy (first term), interaction energy 
(second) and thermal energy has a direct consequence for the phase boundary Tcr{n) of 
Coulomb crystals, cf. Fig. HI While for a classical crystal, the slope of the boundary is 
always positive, dTcr{n)/dn > 0, given by a constant value of F, for quantum crystals, 
there exists a maximum value of the temperature, T^°'^, where the slope changes sign. 
For densities to the left of the maximum the phase boundary is dominated by "normal" , i.e. 
thermal melting, whereas for densities exceeding the value of the maximum, by a competition 
of quantum kinetic and interaction energy. For sufficiently large densities (with decreasing 
Vs) quantum melting is observed, even at zero temperature, cf. Fig. HI The corresponding 
critical values of the Brueckner parameter of a Coulomb OCP at T = are rf ^ 100(160) 
in 3D and rf ^ 37 in 2D for fermions (bosons) 2l|, |45|] and references therein. These values 
are still under investigation. Also, generalization of the results to a quantum Yukawa OCP 



has only recently been attempted, see Ref. 



461] and refences therein. 
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FIG. 4: Left: Melting temperature of small 3D spherical Yukawa crystals versus particle number. 
Right: Phase diagram of the mesoscopic 2D Wigner crystal for different particle numbers A^. OM 
(RM) denotes the boundary of orientational (radial) melting. Here the dimensionless density n 
and temperature T are defined as n = \/2lQ/rQ = (aB/ro)"*^/^ ~ Vg ^^'^ and T = kn T/En, where 



/q = h/{mujQ), Eq = e /ebrQ with tq given by e /e^ro = muj r^/l, from Refs. 



21 



43|. 



Finite trapped quantum plasmas show the same general behavior as a macroscopic quan- 
tum OCP and, in addition, finite size effects as in case of the classical crystals in traps. 
As a consequence, the crystal phase boundary is strongly A^— dependent, as can be seen for 
the 2D case in Fig. HI Further, in 2D the competition of hexagonal (bulk) symmetry and 
spherical symmetry induced by the trap leads to possible additional phases, both in classical 
and quantum trapped plasmas. The most prominent one is a partially ordered phase where 
particle ordering occurs within each shell, but no order of different shells with respect to each 
other exists. Only at significantly larger values of the coupling parameter the orientational 
fluctuations freeze out (orientational freezing or melting, "OM" ) and the crystal enters the 
fully ordered phase, cf. Fig. HI The location of this phase boundary is strongly dependent 
on the crystal symmetry and may vary with by many orders of magnitude |121[]. In 3D 
trapped plasmas no radially melted phase is observed because there is generally a much 
larger energy barrier for intershell rotations. 



B. Critical values of the distance fluctuations 

The appearance of different coupling parameters in the case of classical and quantum 
plasmas makes it very difficult to construct a joint phase diagram of Coulomb crystals. 
An alternative approach to the crystal phase boundary uses, as the starting point, the 
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magnitude of the relative interparticle distance fluctuations of the particles around their 
lattice positions. Expanding, as in Section IV Al the total energy fluctuations AE in a Taylor 
series up to second order and diagonalizing the result allows to express AE as a superposition 
of d ■ {N — 2) relative normal modes. For this system of independent ID quantum harmonic 
oscillators with the phonon modes uj\{q) of polarization A all thermodynamic properties at 
a given temperature T are known. For example, the thermodynamic average of the distance 
fluctuations {6x'^) = (x^) — (a;)^ is given by |2| 

For a macroscopic classical OCP, fx{q,T) 'ntf^^ average over the phonon spec- 

trum yields, in case of a bcc crystal, {Sx"^) = 12.973 f^/F. The result for the relative distance 
fluctuations Urei = {Sx"^) /r^ normalized to the nearest neighbor distance, tq = (37r^)^/^f, 
is 

/ 12.973 1 ^^^^ _ 
<i. = y(^f — 0-155, (7) 

where the last number is the critical value obtained by using F = 175. 

Analogously, we obtain for a quantum OCP bcc crystal at zero temperature, where 

/A(g,r)^i, 

(5^^) = ^«-ir.?aL, (8) 
with «_! = (^) denoting the moment of order minus one of the phonon spectrum which 
equals 2.7986 for a bcc crystal [47|]. This yields, for the relative distance fluctuations 



<^ = \/^^ 0-28 (0-249), (9) 

where the last number is the critical value for fermions (bosons), using r'^^ = 100 (160). 
Note that these fluctuations are mainly due to quantum diffraction effects, i.e. the flnite 
extension of the particle wave functions. Spin effects (quantum exchange) play a minor role 
for the location of the crystal phase boundary which is clear since, in the crystal state, the 
wave function overlap has to be small. Nevertheless, the physical properties of crystals of 
bosons maybe essentially different from the one of fermions. The reason is that interacting 
bosons may show superfluid behavior which may even persist in the crystal phase. This 



state is called a supersolid and was predicted thirty years ago 48|, |49|, |50| and was recently 



observed in PIMC simulations of trapped bosonic plasma crystals |51|- 
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Equations ([7]) and (Q are very useful as they establish the relation between relative 
distance fluctuations and the relevant coupling parameter in the two limiting cases of classical 
and quantum plasmas. To connect the two limits along the whole phase boundary, cf. Fig. 
m one has to use the full phonon spectrum, Eq. (jH]), without expansion of the function fx- 
The temperature and density dependence of Urei remains an open question although some 
interpolations have been attempted, see e.g. Ref. j2|. Further improvements, in particular, 



in the quantum regime, may require to include anharmonic corrections, e.g. Ref. [52], since 
their finite extension lets the particles explore ranges of the potential energy which cannot 
be approximated by a parabola, e.g. Ref. [46 1. 



VI. UNCONFINED TWO-COMPONENT COULOMB CRYSTALS 

As discussed in the introduction, crystal formation in two-component plasmas (TCP) 
competes with bound state formation. One may, therefore, ask whether there exist parame- 
ters where CC exist and, at the same time. Coulomb bound states are ionized. In comparison 
to an OCP, in a TCP, we have at our disposal two additional parameters to realize these two 
conditions: the mass ratio M = mh/nie and charge ratio Z = Ch/ee (in a non-equilibrium 
mass-asymmetric plasma there is further the possibility of different temperatures of the com- 
ponents ^5i|). The first requirement is obvious: the heavy component (ions or holes) has to 
be sufficiently strongly correlated such that it can form an OCP Wigner crystal. The second 
condition is that electrons have sufficiently high kinetic energy to escape the ionic binding 
potential. For classical electrons this requires a sufficiently high temperature whereas in a 
quantum plasma ionization is possible when electron wave functions of neighboring atoms 
start to overlap - this leads to tunnel ionization (Mott effect) which occurs at a sufficiently 
high density. In summary, we find two alternative sets of conditions 

Fj > F"' and ^fc^Te > E^, classical case, (10) 
'Tsi > "rf and r^e < quantum case, (11) 

where in 3D r^°" ^ 1.2. The phase boundary of the Coulomb crystal of the heavy particles 
can be obtained using the harmonic lattice theory results of Sec. IV B[ For the quantum 
case, we may use, for {SxD, Eq. (jS]) and express the nearest neighbor distance of the heavy 
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particles, r^h, by the one of the electrons 



2 _ \"-^h,/ _ " ' sh "'Bh (10) 



2a^ Z^^rlal' 



where atcc = (Svr^)^/^. Assuming that, at the phase boundary, the critical value of the 
fluctuations is given by the OCP result, Eq. (ED, and r,e = r^°" we readily obtain the 
existence conditions of a CC of fermionic (bosonic) ions in a two- component plasma: 



j^^cr^4/3 ^ 3 (132.8). This agrees with the result of Ref. |53j where it was obtained 



from a different derivation. Thus crystallization requires a minimum mass ratio M between 
heavy and light particles. This condition is fulfilled for compact dwarf stars where a crys- 
tal of carbon and oxygen nuclei (fully ionized atoms) is expected to exist 2|. Further 
candidates are crystals of protons which was recently confirmed by PIMC simulations, cf. 
Refs. 



mm, 

particles, see Ref. [53j. Both systems might be accessible in laboratory 
experiments in the near future. Another area where such two-component CC should be ob 



servable are electron- hole plasmas in intermediate valence semiconductors, see Refs. 53l. |55| 



where one could also verify the critical value of M experimentally although values of M 
as large 80 exist only in some special materials. Another promising candidate are charge 
asymmetric bilayers where hole crystallization is expected to occur already for M < 10 [56] 
which is due to the 2D confinement of the particles. 

The analytical predictions of heavy particle crystallization in a TCP have been verified 
by PIMC simulations where both electrons and heavy particles have been treated fully 



quantum mechanically [53|, |55|, |57|] . As can be seen in Fig. [5l with increasing M, indeed hole 
localization becomes more pronounced and, between M = 50 and M = 100, a transition 
to crystal-like behavior is observed. A quantitative analysis based on the relative distance 
fluctuations of the heavy particles, left part of Fig. Ej confirms that the liquid-solid transition 
takes place around M ~ 80. This is a novel kind of quantum phase transition, where melting 
occurs at constant temperature and density - by "changing" the heavy particle mass. The 
phase diagram of the two-component CC is sketched in Fig.[6]for the two values M = 100 and 
M = 200. The larger M the more extended is the crystal phase in the density-temperature 
plane. The crystal phase is bounded from above by the (green dashed) line Th = P^ and from 
the right (high densities) by the (vertical green dashed) line Vgh = rf. This is the simplest 
approximation where the influence of the electrons on the heavy particle interaction has 



been neglected. Improvements require the inclusion of screening effects 
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46], as discussed 



FIG. 5: (Color) Snapshots of a Coulomb crystal of heavy particles (red clouds) embedded into a 
Fermi gas of electrons (yellow) in a macroscopic two-component (neutral) plasma (spin averaged 
results) for mass ratio M = 12 (top left), M = 50 (top right), M = 100 (bottom left), M = 400 
(bottom right). The density corresponds to r^g = 0.64, the temperature is = Th = 0.06Er. 
First-principle two-component PIMC simulations. 

above, this leads to a destabilization of the crystal. At the same time, the heavy particle 
crystal also influences the spatial distribution of the electrons which stabilizes the crystal 
compared to the OCP case. Thus, there exist two competing effects for the crystal stability. 
A detailed comparison of the crystal phase diagram in an OCP and a TCP, therefore, remains 



an interesting still open question. Finally, it has been predicted by Abrikosov [58|] that, in 
the presence of a hole crystal, the electrons should tend to form Cooper pairs, i.e. exhibit 
superconductivity which yet remains to be verified experimentally. 



VII. CONCLUSIONS 



In this paper we have given an overview on strong correlation effects in classical and 
quantum plasmas, in particular on Coulomb (Wigner) crystallization. We have discussed 

15 



0.25 





- 




" --M= 12 






M= 25 




0.20 


\^ ~' "m= 50 






\ M = 75 




0.15 


\ M = 100 






\ /m = 150 




0.10 










• • ■ 


0.05 




• 



100 

Mass ratio M 







1 K 




Electron- i 


Classical hole liquid !• 




hole 1 








bound \ 












Hole crystal 


Quantum 




C 


M=200 


hole liquid 



1/rse 



FIG. 6: (Color) Left: mean-square relative heavy particle distance fluctuations versus mass ratio M 
for Te = 0.096 and Tgf. — 0.63. Symbols are simulation results, the line is the best fit [55]. Right: 
Qualitative phase diagram of a Coulomb crystal of heavy particles ("holes" ) in a macroscopic 



two-component (neutral) plasma. Tg = ^ksT / Eji. Taken from Refs. [5i 



1,Q. 



the possible occurences of Coulomb crystals, first, in trapped one-component plasmas and, 
second, in two-component neutral plasmas. The conditions for crystal formation have been 
summarized in terms of known critical values for the coupling parameters as well as in terms 
of critical values of the relative interparticle distance fluctuations. Using the data for the 
critical parameters it is possible to construct the phase diagram of strongly coupled Coulomb 
matter which was discussed for two cases: mesoscopic classical and quantum plasmas in a 
parabolic 2D trap and two-component mass-asymmetric plasmas. 
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